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SUMMARY 


Aerodynamic equations with unsteady effects were formulated for an 
aircraft in one-degree-of-freedom, small-amplitude, harmonic motion. 
These equations were used as a model for aerodynamic parameter 
estimation from wind tunnel oscillatory data. The estimation algorithm 
was based on nonlinear least squares and was applied in three examples to 
the oscillatory data in pitch and roll of a 70° triangular wing and an X-31 
model, and in-sideslip oscillatory data of the High Incidence Research 
Model 2 (HIRM 2). All three examples indicated that a model using a 
simple indicial function can explain unsteady effects observed in measured 
data. The accuracy of the estimated parameters and model verification 
were strongly influenced by the number of data points with respect to the 
number of unknown parameters. 
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parameter in indicial function 
wing span, m 

parameter in indicial function, 1/sec 
rolling- and pitching-moment coefficient 
normal-force coefficient 
indicial function 


covariance matrix of estimated parameters 


parameter in indicial function 
wing mean aerodynamic chord, m 
deficiency function 

cost function 
reduced frequency, k = 

c b 

characteristic length, l = — ox l- — , m 
& 2 2 


matrix in normal equations 
number of frequencies 
number of angles of attack 
roll rate and pitch rate, rad/sec 
wing area, m 2 
variance 

a. ) standard error 

b. ) parameter in Laplace transform (Appendix A) 
time constant, sec 

time, sec 

aerodynamic derivatives in model for in-phase and out-of- 
phase component of oscillatory data, respectively 
in-phase and out-of-phase components of oscillatory data, 
respectively 


term in model equation (6), 


z u = 


1+T^ 2 


term in model equation (7), 


z„ = 


II 


1+ t 2 & 2 


2 



a angle of attack, rad or deg 

P sideslip angle, rad or deg 

S c canard deflection, deg 

0 vector of unknown parameters 

t time delay, sec 

Tj nondimensional time constant, = — 

b]£ 

</) roll angle, deg 

co angular frequency, rad/sec 

Superscript: 

estimated value 


Superscript over aerodynamic derivative: 

oscillatory data in pitch or roll 
oscillatory data in heave or sideslip 


Subscript: 

A amplitude 

0 nominal value 


Matrix exponent: 

-1 inverse matrix 

T transpose matrix 


Aerodynamic derivatives 

dC A 


C, ( oo) = . 


da 


for A = l,m, or n and a = 


qc 2 a fib 
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INTRODUCTION 


Reference 1 presents a short theoretical study of aircraft 
aerodynamic model equations with unsteady effects. In one of the 
examples in this study, a formulation of unsteady aerodynamics is applied 
to small-amplitude, one-degree-of-freedom harmonic motion of an aircraft 
about one of its body axes. The mathematical models developed can be 
directly used in the analysis of wind tunnel data obtained from forced 
oscillatory tests. During these tests, a model is forced to oscillate in the 
tunnel airstream about any single model body axis at a specified angular 
frequency and amplitude. From the measured aerodynamic forces and 
moments, the in-phase and out-of-phase components of the aerodynamic 
coefficients are obtained (see reference 2). It has been a long-standing 
practice to analyze the oscillatory data by formulating the aerodynamic 
model equations with constant parameters known as stability derivatives. 

In this formulation, the in-phase component is a combination of a static 
derivative (such as C ma ) and a rotational acceleration derivative (such as 

<v whereas the out-of-phase component is a combination of a purely- 

rotary derivative (such as C m ^ ) and a translation acceleration derivative 

(such as C m& ). 

The results from forced-oscillatory tests with a given amplitude and 
at a given angle of attack show that the resulting combinations of 
derivatives very often depend on the frequency of the oscillations. This 
dependence contradicts the basic assumption about the time-invariance of 
the stability derivatives. The effect of frequency on the aerodynamic 
parameters is explained by a proposed formulation of the aerodynamic 
model equations with unsteady terms. In addition, this formulation 
separates the sums of derivatives, mentioned above, into two terms, one 
representing either the static or purely-rotary derivatives, and the other, 
the unsteady effects. 

Attempts to analyze the wind tunnel oscillatory data with the 
inclusion of unsteady aerodynamic effects are reported in references 4 and 
5. The purpose of this report is to develop a method for estimation of 
unknown parameters in mathematical models postulated for wind tunnel 
data from small amplitude oscillatory testing and demonstrate this method 
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in three examples. In the first example, the forced oscillatory data in pitch 
and roll of a triangular wing are analyzed. The second example uses the 
oscillatory data in pitch and roll of a model of the X-31 experimental 
aircraft. The third example deals with acceleration-in-sideslip data of a 
high incidence research model (HIRM). In this test, the model oscillates in 
sideslip at a specified angular velocity and amplitude (see reference 5). 

From the measured data, the aerodynamic derivatives corresponding to 
sideslip angle and its rate are calculated. 


PARAMETER ESTIMATION 


A parameter estimation algorithm is developed for the analysis of the 
normal force obtained from wind tunnel oscillatory data in pitch. 
Modification of this algorithm for oscillatory data in roll and yaw, and 
oscillatory data in heave and sideslip should be straightforward. As 
pointed out in reference 1, the equation for the normal force can be 
formulated as 


C N (t) = C N (oo)a{t)- \F a {t-T)-^-a{T)dT + ^-C N (°°)q(t) (1) 

J dr V Q 

0 

where 

Qv a (°°) is the rate of change of C N with a in steady flow evaluated at 

q = 0. This term corresponds to the stability derivative — — . 

da 

^N q (°°) is the rate of change of Cjy with q in steady flow evaluated at a 


equal to its mean value during the oscillatory motion. This term 

dC 

corresponds to the stability derivative — 

n Q c 


F a (t) is the deficiency function which is defined as the difference 
between Cjiy a (°o)and the indicial function Its analytical 

form, proposed in reference 1, is 

F a (t) = ae~ blt 

which follows from the postulated indicial function 
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( 2 ) 


Cn k (£) = a| 1 - e~ blt j + c 
= C iV a ( 00 )-- p, a^) 

The steady-state form of equation (1) for harmonic changes in a(t) is given 
as 


where 


— I — 

Cjy (t ) = C]$ a a & sin(ft# ) + — Cj^ a^co cos (ait) 


C N a and are the in-phase and out-of-phase components, 
respectively, 

a & is the amplitude of the oscillations, 
ft) is the angular frequency, 
l is the characteristic length, and 
V is the airspeed. 


(3) 


As follows from the development in reference 1, the mathematical models 
for the components of Cpj(t) have the form 


C*,=C„ 9 M-a I -^ r 


(4) 

(5) 


where is related to the parameter as 


V 



and k is the reduced frequency 


k 


(M 

V 


The second term on the right side of equations (4) and (5) represent the 
unsteady counterparts of and -Cjy\ j respectively. The parameter a 

accounts for the variation of the unsteady terms with the angle of attack. 
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From the experiment, the in-phase and out-of-phase components are 
usually obtained for different values of the angle of attack and reduced 
frequency while keeping the amplitude of the oscillations constant, that is 


Qv a { a i’kj) 
Qv { a i’kj) 


i = 1 > 2 , 

7 = 1 , 2 , 


n 

m 


Then equations (4) and (5) can be expressed as 

u ji ~ u i f ( ) — a i z uj f{®i ) (6) 

Vji = v i ~ ®i z Vjf(^i ) (7) 

where 


u i = c N a (°°; ), = C Nq (oo; ), 

x 2 k 2 

1 7 _ T 1 

ZU -> = l + X 2 k 2 ’ Zv j = l + T 2 k 2 ’ 

/(«i) = 1 

The definition of z u . and z v . for oscillatory data in roll and yaw are 

identical to those for the data in pitch, but for the rolling oscillations, 
f(a) = sin (a) and for the yawing oscillations, f(a) = ±cos(a) where the 
minus sign applies to the out-of-phase components. 

In equations (6) and (7) there are, in general, 3n + 1 unknown 
parameters: Mj, v t , a i} and t*. They can be estimated from experimental 
data uji and uji by minimizing the cost function 



This cost function is nonlinear in the parameter T\. To formulate a linear 
estimation problem, equations (6) and (7) are linearized about some 
nominal values of all unknown parameters. The linearized cost function 
takes the form 
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where 



and the index 0 indicates nominal values. 

Following the minimization of (9), the normal equations for unknown 
parameters can be formulated as 


where 


Ad = -M 


_i 9J 

90 


0 = 0 o 


( 10 ) 


AO = [A zi; A V} Adi Ati]" 1 


and M is the matrix formed by partial derivatives 


9u 


Ji 


90 


and 


9v 


Ji 


0 = 0 n 


90 


0=0n 


Parameter estimation using equation (10) is an iterative process where the 
estimates are obtained as 


0 r + 1 = 0 r + ( 11 ) 

and where the index r indicates the r iteration. The parameter 
covariance matrix is estimated as 

Cov(0) = s 2 M _1 (12) 


with the variance estimate 



(13) 



2nm-(3n + 1) 


For the estimation of all parameters in (6) and (7), the number of 
measured data must be equal to or greater than 3n + 1 . This means that 
the experiment must be repeated for at least two different frequencies, that 
is m> 2. The estimation algorithm can be modified by the inclusion of 
C N J oo * aj) from static wind tunnel tests. Then the measured data uji will 

be replaced by [uji -u^j. The number of unknown parameters will be 2n + 1 

and again, the condition m > 2 will have to be met. If, on the other hand, 
only the out-of-phase components vji are available, then the number of 

data points is nm, the number of unknown parameters is 2n + 1, and the 
condition for the number of frequencies is m > 3 . 

The model for the normal force of an aircraft performing one-degree- 
of-freedom harmonic oscillations in heave is formulated in Appendix A as 

i ~ 

C N (t) = C^ a aAsm(ojt) + — Cj^ (k a J ^o)cos((ot) (14) 

where 


t 2 & 2 


CN n = ~ a 


1+ t 2 & 2 


(15) 

(16) 


Following the previous development, equations (15) and (16) can be 
expressed in the form 


“jfi = u i~ a i z uj 
Vji ~ ~ a i z vj 


i = 1,2,.. 

7 = 1 * 2 ... 


,n 


C17)— 

(18) 


The algorithm for estimation of unknown parameters Uj, a j, and Tj can be 
obtained in the same way as for the parameters in equations (6) and (7). 

The estimation algorithm using oscillatory data in roll and sideslip 
follows from the model developed in reference 1 
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T 2 * 2 


Q =Q (oo)sm(a) a 1 9 9 sin(a) 
p p 1+rfk 2 

(19) 

Cl p =Cl p (°° ) a 1+T 2 A 2 Sinta) 

(20) 


and from the model developed in Appendix A 




(oo)-a 


rfe 2 

1+ r 2 k 2 


(21) 



-a 


*1 

1+ r 2 k 2 


( 22 ) 


Similar to the expressions for the oscillatory data in pitch, the second term 
on the right side of equations (19) and (20) are the unsteady counterparts of 

Q 

k Ci p and — Q sin(a), respectively. 


EXAMPLES 

The estimation algorithm developed is applied to measured 
oscillatory data in three examples. In the first example, parameters of a 
triangular wing subjected to forced oscillations in pitch and roll are 
estimated. The second example again presents the analysis of oscillatory 
data in pitch and roll this time, however, measured for a 19-percent-scale 
model of the X-31 experimental aircraft. Finally, in the third example, 
parameter estimates are obtained from data due to acceleration in sideslip 
of the High Incidence Research Model 2 (HIRM 2). A sketch of the last two 
models is given in figure 1 together with the numerical values of their basic 
geometric characteristics. 

Forced Oscillations of 70° Triangular Wing 

The oscillatory data of a triangular wing were obtained from 
reference 5. The model had the span, b = 0.90 m, mean aerodynamic 
chord, c = 0.824 m, and sweep angle of the leading edge of 70°. The model 
was tested in pitch and roll at two center-of-gravity (axis of rotation) 
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locations, 0.25 c and 0.50 c, and at angles of attack between -4° and 57°. The 
frequencies varied between 0.5 and 2.5 Hz and the amplitudes of the 
oscillations were 3° and 5°. The experiment was conducted at the speed of 
45 m/sec corresponding to a mean-chord Reynolds number of about 2x10®. 
Reference 5 contains only a limited number of test data: the in-phase and 
out-of-phase components of the normal force and the out-of-phase 
components of the pitching and rolling moment. For parameter 
estimation, the data with the amplitude of 3° and c.g. location of 0.50 c were 
chosen. The following three sets of data were used in the analysis: 

1. , kj ) and Cjy at three frequencies and eight values of the 

angle of attack between 27° and 56 °, 

2. C m ^ ( oci,kj ) at five frequencies and the same values of angle of attack as 
in the preceding set, 

3. C/ (a; , kj ) at five frequencies and fourteen values of the angle of attack 
between 5° and 57°. 

The measured data are plotted against the angle of attack with the reduced 
frequency as a parameter in figures 2, 4, 5, 7, and 8. 

The parameter estimation algorithm for data set C ^ a , was based 

on model equations (4) and (5); for data of C m<? , on equation (5); and for data 
of Q , on equation (20). The estimates of the nondimensional time delay, 

V 1 

Ti, and computed values of parameters 6^ = and Ti = — for all data 

b x 

sets are given in table I. The values of Ti are similar for the normal-force 
and pitching-moment indicial functions, and only slightly higher for the 
rolling-moment indicial function. The increased number of data points due 
to increased number of frequencies from 3 to 5 improved substantially the 
accuracy of the estimated parameter Tp 

Estimates of parameter a in the indicial function and derivatives 
C 'N a (°°), (°°), (°°), and C[ p (°o) are plotted in figures 3, 6, and 9. 

The minimum and maximum values of their standard errors are 
presented in table II. All the estimates mentioned are obtained from data 
at three frequencies. Contrary to the estimates of Tj, increasing the 
number of frequencies from 3 to 5 had only a small effect on the estimates of 
the remaining parameters and their covariances. 
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In figure 3, the estimates of Cjv a (°°) are compared with the results of 

a static wind tunnel test. The agreement between estimates and measured 
data is very good. This agreement can be interpreted as a confidence in the 
postulated model for C^ a and • The model was further verified by a 

comparison of the measured and estimated values of the oscillatory data 
and by checking the model prediction capabilities. In figures 2, 4, and 7 the 
measured data are plotted together with those estimated from model 
equations after substituting the parameter estimates. The agreement 
between measurement and estimation is very good. In figure 5, the 
measured values of C m are compared with those predicted at two 

frequencies not used in the parameter estimation. As in the previous 
cases, the predictions agree very well with the measured data. Similar 
agreement can be seen in figure 8 where the measured and predicted 
values of Cj are shown. 

l p 

If the measured data are obtained only for a limited number of 
frequencies, good agreement between measured and estimated points in 
itself cannot be considered as a verification of the model adequacy. As 
follows from equations (4) and (5), at each selected angle of attack there are 
four parameters which are to be estimated from 2m data points of the in- 
phase and out-of-phase components. If m- 2, only the constraints a.) Tj 
has the same value for all a { - and b.) aj has the same values for C^ a and 

Ctfq prevented complete agreement between measured and estimated data. 

Accurate parameter estimation and model verification would, therefore, 
require measurement of the in-phase and out-of-phase components at an 
increased number of frequencies (four or more for estimation, one for 
verification) and static wind tunnel testing for obtaining aerodynamic 
derivatives in the model equations. 

Forced Oscillations of X-31 Model 

The forced-oscillation test on the model of the X-31 aircraft were 
conducted at the 30-by-60 foot wind tunnel at NASA Langley Research 
Center. During the test, the dynamic pressure at the tunnel was 10 lb/ft^ 
(478.8 Pa) and the Reynolds number referred to c was 1.37x10®. From wind 
tunnel results, three sets of data were selected for parameter estimation: 
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1. C/y a and at two frequencies of 0.8 and 1.2 Hz and 0° < a < 87.5°, 

2. C ma and C m ^ at the same frequencies and angles of attack, 

3. and at two frequencies of 0.4 and 0.6 Hz and 25° < a < 50°. 

The amplitude of oscillations was 5°. The positions of the control surfaces 
were as follows: the trailing edge, 0°; leading edge, 40° (inboard) and 32° 
(outboard); and canard, -40°. Some of the measured data are available in 
reference 7. 

Similar to the previous example, model equations (4), (5), (19), and (20) 
were used. The measured and estimated values of C^ a and are 

shown in figure 10. The estimated points are almost identical or close to the 
measured data. This closeness is the result of the small number of data 
points for a given angle of attack, as pointed out in the previous example. 

The estimated parameters T\ and computed parameters b\ and 7\ are 

given in table II together with the same parameters estimated from the 
remaining two sets of data. The parameters Cjy a (° °), Cjy^ (°°), and a are 

plotted in figure 11. Substantial differences in the derivative Cjy a (<*>) 

estimated from the oscillatory and static measurement and large scatter in 
the estimates of Cpj (°°) around the fitted curve are apparent. 

Estimated parameters from data sets C ma , C m ^ and , and Q ■ are 

presented in figures 12 to 15. These estimates have problems similar to 
those mentioned before. In addition, the review of Cj at various 

L p 

frequencies, amplitudes, and configurations presented in figure 16 leads to 
the conclusion that the estimates of Q (°°) in figure 15, for angles of attack 

between 25° and 35°, have values far from those expected. Low parameter 
accuracy and differences in derivatives C/y a (°°), (°°), and (°°) from 

oscillatory and static measurements could be caused by a small number of 
data points, measured data accuracy, and by modeling errors in postulated — 
models used in parameter estimation. For obtaining more accurate 
results, it would be necessary to have the oscillatory data (both in-phase and 
out-of-phase) at more frequencies and more values of the angle of attack 
within the selected range. 
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In-sideslip Oscillations of HIRM 2 

A description of the HIRM 2, testing method and conditions of the 
test, and some of the results are given in reference 5. The experiment was 
conducted in the 13-by-19 foot low speed wind tunnel at the Defense 
Research Agency (formerly Royal Aircraft Establishment) in Bedford, 
United Kingdom. The wind speed was approximately 30 m/sec which 
corresponds to a Reynolds number of about 1.1x10®. For parameter 
estimation, two sets of data were selected: 

1. C/^ and Q at three frequencies of 2, 3, and 4 Hz and sixteen angles of 

attack between 0° and 42°, canard off, 

2. and C^. at the same test conditions, but with the canard on at 0° 

setting. 

The model for the measured data is represented by equations (21) and 
(22) with three unknown parameters (°°), a, and X\. The estimated and 

computed parameters Tj., &]_, and T\ are summarized in table IV. These 
parameters have the same values for both configurations. The time 
constant is about one-third of that in the two previous examples. This 
means that the unsteady effect in sideslip oscillations of the HIRM 2 is less 
pronounced that in the oscillatory motions of the triangular wing and the 
X-31 model. The estimated values of and C/ presented in figures 17, 

18, 20, and 21 are mostly close to the measured data which might be 
considered as the first indication of an adequate model for the given data. 

The remaining two parameters are plotted in figure 19 and 22. Rather 
large differences exist between estimates of Q (°°) from oscillatory and 

static wind tunnel data obtained from reference 8. The reason for this 
disagreement was not found. No further model verification was possible 
due to the limited amount of measured data points. Plots of the parameter 
a indicate that, for the angles of attack between 0° and 25°, the effect of 
unsteady aerodynamics is very small. Therefore, for better understanding 
of the unsteady phenomenon and its modeling, it would be necessary to 
increase the number of data points by increasing the number of frequencies 
and selected values of angle of attack within the range from 25° to 50°. 
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CONCLUDING REMARKS 


Using nonlinear least squares, an estimation algorithm for aircraft 
aerodynamic parameter estimation from wind tunnel oscillatory data was 
developed. Models of an aircraft in one-degree-of-freedom, small- 
amplitude, harmonic motion included unsteady terms in the form of 
indicial functions. In formulating analytical form of indicial functions for 
this type of data analysis two conflicting requirements must be addressed: 
parameter estimation requires a simple model with a small number of 
parameters in order to insure their reliability; on the other hand, a simple 
model will not completely explain the rich and complex phenomena at 
various scales associated with unsteady and separated flow during 
oscillatory or transient motion. The indicial functions were postulated as 
simple exponential forms where the unknown parameters included 
aerodynamic derivative, the exponent and multiplication term. It is 
important to realize that the model proposed in this work should be used in 
the analysis of experimental oscillatory data where the effect of frequency at 
a given nominal angle of attack and Reynolds number is considered. For 
different applications the proposed form of indicial functions should be 
carefully considered. 

The estimation procedure was applied in three examples to the 
oscillatory data in pitch and roll of a 70° triangular wing and an X-31 
model, and to in-sideslip oscillatory data of the High Incidence Research 
Model 2 (HIRM 2). From postulated models and examples, the following 
conclusions can be drawn: 

1. An aerodynamic model with a simple form of the indicial function 
included can "explain" unsteady effects observed in small-amplitude, 
oscillatory wind tunnel data; 

2. the accuracy of estimated parameters and model verification could be 
strongly influenced by the number of measured data points and their 
accuracy; 

3. future wind tunnel experiments intended for parameter estimation 
should include oscillatory data (both in-phase and out-of-phase 
components) at a large number of frequencies, five or more, and at an 
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increased number of angles of attack in the region where unsteady 
effects can be expected; 

4. the wind tunnel experiment should also include measurements for 
obtaining static stability derivatives and purely-rotary derivatives for a 
comparison with their estimates from oscillatory data. 
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APPENDIX A 


HARMONIC OSCILLATORY MOTION IN HEAVE AND SIDESLIP 

Using equations (1) and (2), the expression for the normal force of an 
aircraft performing a one-degree-of-freedom (o.d.f.) oscillatory motion in 
heave has the form 


t . 

C^(t) = Cjy (oo)a(t)-a fe -6l (* -T ) — a(x)dx 
a J dr 


(A.l) 


Applying the Laplace transform to (A.l), the expression for the normal- 
force coefficient is obtained as 


Ctf(s) = 


f r as A 


a(s ) 




where, for simplicity, Qy a = (°o). 

Using a complex expression for harmonic changes in a(t), that is 

a(t) = oc A e l03t = «A[cos(ftrf) + isin(G#)] 
and replacing s by io), the steady-state solution to equation (A.2) is 


(A.2) 


Ctf(t) = 


C N„ ~a 


a r 


a b 2 + co 2 


a A sin{cot)-a 2 1 ^ a A (ocos(ajt) (A.3) 


&i + CO 


The introduction of reduced frequency k = ~^- and nondimensional time 

V 

constant x^ = — yields 
b\( 


where 


Cjyit) = C^ a a A sin(cot) + C^ di a A kcos((jot) 


(A. 4) 


c N a = C Na (oo)- a —l 


xH 2 


l+x 2 k 2 


(A. 5) 


&Na = ~a 


1+ x 2 k 2 


(A.6) 
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Similar expressions can be obtained for the pitching moment coefficient 

Based on the analogy with the preceding case, the rolling-moment 
coefficient of an aircraft performing a o.d.f. oscillatory motion in sideslip 
has the form 


Ci(t) = C h P A sin(cut) + Ci^f3 A kcos(cot) (A.7) 

where 




(°° )~a 


T ^ 2 

1+T 2 ^ 2 


(A.8) 



-a 


1+ T 2 ^ 2 


The expressions for the side-force and yawing-moment coefficient are 
identical to those in equations (A.7) through (A.9). 


(A.9) 
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Table I. Estimated and computed parameters. 70° triangular wing. 


Measured 

Number of 

Parameter 

data 

frequencies 

*1 

6 1; sec 1 

7l,sec 

c Na ,c Nq 

3 

24.3 ± 0.29* 

4.5 

0.22 

Cm, 

3 

5 

30. ±6.7 
29. ±3.0 

3.6 

3.8 

0.28 

0.27 

c, 

l p 

3 

40. ± 12. 

2.5 

0.40 

5 

31. ±2.3 

3.2 

0.31 


the second value indicates the standard error 


Table II. Minimum and maximum values of standard errors of 
estimated parameter. 70° triangular wing. 


Measured 

data 

Standard error 

s(d) 


min 

max 

min 

max 

Cn 

iy a 



1 

0.42 


0.040 

0.050 

■ 


N 

ly q 



■ 

0.46 

Cm, 

0.036 

0.091 

0.37 

0.95 

Cl 

L p 

0.008 

0.036 

0.059 

0.20 


■3 1 

for index A a = N a ,N q ,m q , or l p 
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Table III. Estimated and computed parameters. X-31 model. 


Measured 

data 

Parameter 

n 

■119H 

T^sec 

CN a .c Nq 

23.4 ±0.33* 

3.3 

0.30 

C m a >Cm q 

13. ± 1.4 

6.0 

0.17 


10. ± 3.0 

4.2 

0.24 


the second value indicates the standard error 


Table IV. Estimated and computed parameters. HIRM 2. 


Measured 

Parameter 

data 



Tj_,sec 

canard off 

4.6 ±0.38* 

11.8 

0.085 

canard on 

4.5 ±0.41 

12.0 

0.083 


the second value indicates the standard error 


21 


















Figure 1. Two-view sketch of X-31 model and HIRM 2. 
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Figure 2. Measured and estimated in-phase and out-of-phase 
components of normal force. 70° triangular w in g 






Figure 4. Measured and estimated out-of-phase components in pitching 

moment. 70° triangular wing. 



Figure 5. Measured and predicted out-of- phase components of pitching 

moment. 70° triangular wing. 
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Figure 7. Measured and estimated out-of-phase components of rolling 

moment. 70° triangular wing. 



Figure 8. Measured and predicted out-of-phase components of rolling 

moment. 70° triangular wing. 
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a, deg. 

Figure 9. Estimated parameters of rolling-moment component. 

70° triangular wing. 
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Figure 13. 
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Estimated parameters of pitching-moment components. 
X-31 model. 


32 







a, deg. 

Figure 14. Measured and estimated in-phase and out-of-phase 
components of pitching moment. X-31 model. 







Figure 16. Effect of canard setting, frequency, and amplitude of 
oscillations on out-of-phase component of rolling moment. X-31 model 
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